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Anomalous Threshold as the Pivot of Feynman Amplitudes 
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Reduction techniques, Landau singularities and differential equations for Feynman amplitudes are briefly re- 
viewed. 
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1. Reduction techniques and factorization 

A modern version of reduction of Feynman in- 
tegrals [[1] tells us that 

|l-loop n-legs Feynman diagrams | = 

J2b v d (p?,...,p?)+--- (i) 



V 



where D is a partition of {l...n} into 4 non- 
empty sets, P® is the sum of momenta in i G V 
and Do a scalar box. In other words, scalar 
one-loop integrals (up to boxes) form a basis. 
Thus, coefficients in the expansion (B-p etc.) are 
uniquely determined, although some reduction 
method can be more efficient than others. How- 
ever, troublesome points where the numerical sta- 
bility of the result is at stake will always be 
there. What to do in these cases? We can change 
(adapt) bases, or avoid bases (expansion). 

We explain our idea via examples; first, we con- 
sider factorization of Feynman amplitudes, the 
Kershaw theorem of Ref. [ [2]: any Feynman 
diagram is particularly simple when evaluated 
around its anomalous threshold. The singular 
part of a scattering amplitude around its lead- 
ing Landau singularity may be written as an al- 
gebraic product of the scattering amplitudes for 
each vertex of the corresponding Landau graph 
times a certain explicitly determined singularity 
factor which depends only on the type of singu- 
larity (triangle graph, box graph, etc.) and on 
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the masses and spins of the internal particles. 

Let us illustrate the consequences of factoriza- 
tion with one example: define a scalar one-loop 
N -leg integral in n -dimensions as 



A, 



with A„ = (J, 4 n /(in 2 ) and Pi = q -+ 
(po =0). In parametric space we have 



Nn = — 



2\ N-n/2 



r(7v--J K 



v N 



i=i 



X 1 H K 



dxj V f 



2K{ 



i/2-N 



(2) 



X N — A , r U Ai 



Standard notation for A = 1,2 ... is Nq = 
A ,Bo . . .; the superscript n will be dropped un- 
less strictly needed. In order to discuss the pro- 
cedure it is helpful to introduce the following 
quantities: the BST factor [ |3j, B N = L N — 



Pi ■ Pi 



K f N H N 1 K N , the Gram matrix, 
the Caley matrix [ U] 



I H N K A 
1 Kl L N 



It follows that [ [5j B = C/G for any N, where 
C = detM and G = dctH. Landau singu- 
larities are seen as pinches (we assume that 
masses and invariants e R) when we write V N = 
(x — X N ) 1 H (x — X N ) + B N . This realtion in- 
deed shows that B N = is the origin of the pinch 
on the integration contour at the point of co- 
ordinates x = X N ; therefore, if the conditions, 
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B N — and < X N:N _ 1 < ... < X N>1 < 1, 
are satisfied we will have the leading singularity 
of the diagram (hereafter AT). 

Nowadays, the keyword in any reduction proce- 
dure is to avoid inverse powers of Gram determi- 
nants. A common wisdom, but why? The vanish- 
ing of the Gram determinant is the condition for 
the occurrence of non-Landau singularities, con- 
nected with the distorsion of the integration con- 
tour to infinity; furthermore, for complicated di- 
agrams (see Sect. 10 of Ref. [|6]), there may be 
pinching of Landau (C = 0) and non-Landau sin- 
gularities (G = 0), giving rise to a non-Landau 
singularity whose position depends upon the in- 
ternal masses [\7\. 

Given the above properties, the factorization 
of Kershaw theorem [[2] follows. The beauty of 
being at the anomalous threshold is that scalar 
products are frozen as a consequence of the 
Landau equations and the amplitude factorizes. 
Therefore, the AT looks perfect for boundary con- 
ditions, as long as it is inside the physical region. 
Alternatively we may expand and match residues 
at a given AT [HJ. 

Let us consider standard reduction [[5] as com- 
pared to modern techniques [[5] by taking a box 
diagram with q ■ pi in the numerator: 

3 3 

D-p 1 =^D li p 1 -p i = -^D li H u . (3) 

i=l i=l 

A careful application of the standard method 
gives 

1 



H- „- ^ d j , di 



D. 



(<+l) 



D, 



2KiD 0i 



where D is the scalar triangle obtained by re- 
moving propagator i from the box. Therefore we 
obtain 



D-jn 



1 



1 



• di, 



(4) 



without having to introduce G3. Furthermore, 
the coefficient of the scalar Do in the reduction 
is 1/2 (mo - m? — p\). At the AT of the box we 



must have q 



0, (q+pi) 2 + m\ = 0, etc. 



Therefore the coefficient of Dq is fixed by 



2q-px 



222 
= m -m 1 -p u 



(5) 



which is what a careful application of standard 
reduction gives. Note that one gets the coeffi- 
cient without having to require a physical singu- 
larity. In standard reduction for a N -point func- 
tion each, reducible, scalar product in the numer- 
ator is replaced by a difference of propagators plus 
a K -factor. The latter is what is predicted by fac- 
torization at the anomalous threshold; the proce- 
dure is continued and one finds N — 1 point func- 
tions with reducible and also irreducible scalar 
products; for the latter inverse powers of Gram 
determinants remain. 

It is worth noting that starting from six legs 
factorization must be understood as performed 
at some subLeading Landau singularity of the di- 
agram [110). which is easily achieved by using the 
BST-algorithm [[3]. If the derivation is to hold 
we must further require that the leading Landau 
singularity point does not also lie on the Landau 
curve of its sub-graphs. For illustration, consider 
a box in n -dimensions in a region where -B4 ^ 0. 
BST relations allow us to decompose the box in 
a n + 2 -dimensional box plus four n -dimensional 
triangle, D n — » D n+2 © C n . A second application 
gives D n+A © C n+2 © C n . A box in 8 -dimensions 
as well as a triangle in 6 -dimensions cannot de- 
velop a singularity, threfore the subleading singu- 
larities of the original box are given by the leading 
ones of the four triangles obtained by shrinking 
one of the lines in the box to a point. The coef- 
ficients of the decomposition can be found in [[S] 
and the argument can be generalized to arbitrary 
number of legs. 

To summarize, at least in one point we can 
avoid reduction, all integrals are scalar; however, 
we need to have the AT inside the physical region 
Rphys (support of A ± -propagators in R) Since 
this is a rare event we must have a generalization 
of the factorization theorem: prove that the AT, 
even with invariants ^ i? P hys implies a frozen q. 

If a one-loop, iV-legs scalar diagram is singular 
at x = X N £ R then consider N nfl p^ 1 , 

N 

N n -pi = -^Jf(i)prft 
1=1 

N N 

A^"(i) H u at J2 H U X l = — KiAf n - 
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where N n (i) is the same as the scalar integral 
(JV"(1)) but with one power Xi in the numerator, 
and H X = —K: this leads to generalized factor- 
ization since, at the AT, all scalar products are 
replaced by the solution of (q+. . .+Pi) 2 +mf = 0, 
with i = 0, ... ,7V - 1. 

2. Feynman diagrams aroud AT 

In this section we consider a classification of 
physical ATs: for instance, direct calculation 
shows that, for N = 4, there are 14 branches 
in p-(real) space. In general, this classification 
is much easier when we use the Coleman - Nor- 
ton theorem [ HT]. As a consequence of it, in a 
2^2 process, two unstable particles in the ini- 
tial state are needed. Other simple examples of 
physical AT are represented by a) 7*(Z*) — » bbH 
(for a virtuality s > 4 m 2 and M 2 (if) > 4 ) 
and b) from pentagons arising in the reduction 
of the hexagon in 7*(Z*) — > bbvvH (as suggested 
by A. Denner). 

The expansion of Feynman integrals around 
their AT is easy to derive analytically and only 
requires Mellin-Barnes and sector decomposition 
techniques as explained in Ref. [|5|. Examples of 
leading behavior are: for the vertex Co ~ In B3; 

— 1/2 

for the box Dg ~ B A ; for the pentagon Eq ~ 
B^ 1 and no singularity for the hexagon Fq in 
4 dimensions [ [10] ; e.g. Im Co has a logarithmic 
singularity, Re Co has a discontinuity. Here we 
do not consider infrared/collinear configurations 
where we expect an enhancement of the singular 
behavior (in the residues of IR/coll. poles). 

It is worth noting the non-integrable (scalar) 
pentagon singularity which seems to require the 
introduction of complex masses for unstable in- 
ternal particles For integrable singularities 
we always average over a Breit-Wigner of the in- 
variant mass of unstable external particles. 

3. Differential equations 

An interesting feature of factorization at AT is 
the possibility of introducing a differential equa- 
tion with boundary conditions at the AT where 
the amplitude is directly given in terms of scalar 
functions; what we want is an ODE for the full 



amplitude, with real momenta and one boundary 
condition; this requires to find the right variable. 
The advantages of this procedure are given by a 
total absence of reduction and by the extedibility 
to higher loops. 

It is well-known that non-homogeneous sys- 
tems of ODE [[13] are easy to obtain with IBP- 
techniques [ [14] but the non-homogeneous part 
requires (a lot) of additional work; the natural 
alternative would be to introduce PDE. They are 
notoriously much more difficult to handle even if 
homogeneous (compatible) systems of nth-order 
PDE are easy to derive, a fact that has to do 
with the hypergeometric character of one-loop di- 
agrams. It is enough to use Kershaw expansion 
around pseudo-threshold [ [15] and a generaliza- 
tion of Horn-Birkeland-Ore theory [116]. 

4. Diffeomorphisms 

Let us restrict to ODE. To achieve our goal 
we find it most natural to introduce special dif- 
feomorphisms T of the Feynman diagrams. De- 
fine Pi(z) = T lj (z)p j with J2 Pi — J2 Pi — 
and with Ty (0) = ; next we look for a z = 
z AT £ R where the transformed diagram is sin- 
gular. Furthermore, T is physical if maps -D(O) 
onto a D(z) which is singular at z AT £ R and 
Sij — » Sij(z) £ Phys^, where Sij and Sy are in- 
variants; no restriction on sy is required. T is 
unphysical if maps D(0) onto a D{z) which is sin- 
gular at z AT £ R but — > Sij(z) ^ Phys z ; it 
requires restrictions on the original invariants . 

A general solution of our problem is as follows: 
if 3 a diagram D, a transformation T such that 
D{z) = T(z) D with T(0) = I and D(z AT ) singu- 
lar (z AT £ R) then we map D as follows: 

D —* D (z, z AT ) 
D (z, z AT ) = Ti (z, z AT ) D + T 2 (z, z AT ) D(0) 
Ti{0,z AT )=I, T 2 (0,z AT ) = Q 
T\ (z AT , z AT ) = 0, T 2 (z AT , z AT ) = I. (6) 

It is worth mentioning that, in this way, we can 
write a differential equation for the full amplitude 
instead of one for each master integral with differ- 
ent boundary conditions. The interesting feature 
can be summarized as follows: for a given topol- 
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ogy which is candidate to satisfy Coleman - Nor- 
ton (e.g. crossed box in 2 — > 2) we perform the 
transformation in such a way that the new invari- 
ants indeed satisfy the conditions of the theorem; 
for all parent topologies (e.g. direct boxes) we use 
the general mapping described in Eq.©. 

It is straightforward to see how our approach 
is related to the one of differential equations in 
Mandelstam variables: 



dz 



Oij ({s}) = P H 



dsi d 
dsi ' 



dP 



(7) 



jf i 



where r = {dT/dz)T- 1 . 

As an example for a four-point function we con- 
sider one of the crossed diagrams in gg — > 77 with 
a massive loop. The transformation is 

/l-z z \ 
1 - z z 
z 1 - z 
z 1-* 



T = 



V 
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The transformed invariants are M? = z (1 
and 

S=(l-2z) 2 s, T=(l-2z) 2 i, U = u. 



z) u 



The solution of B 4 
integrand is 



which makes singular the 



2 u 2 z A 



s (4? 



4 m 2 s 



r-1) 

t) (4 m 2 s + ui) 



1/2 



The effect of the transformation is simple, we 
have mapped the original box onto a box which 
satisfies the condition stated in Coleman - Norton 
theorem. 

As an example of ODE in z we consider the 
scalar box after the transformation Pi. 4 = p\_ 4 + 

Pi), 



z(pi +P2) and P 2 ,3 = 
£>0 (M) = A„ / d n q 



P2,3 - Z (pi 
1 



n 



i=0,3 



D?(i)=U?(l 
^0 = ^0 (1 

d 



1) 



n« 



2zs 



■ ,2, . 

(2)+P "(4) 



triangles (8) 



Using IBP-techniques (and dropping the super- 
script n) we derive 

D (i) = dj, det R 4 (z AT ) = (9) 



where di contains D or triangles. Introducing 
r — z 2 — z, we obtain 

^D (r) = C^(r) [x(r)D (r)+D rest (r)\ (10) 

where C4 is the corresponding Caley determinant. 
Furthermore, we have 



-Ca 
dr 



-2X(r), 

which leads to the expected solution, 



D (r) 



P s 



(r- 



,1/2 



D les (r) 



(11) 



(12) 



Before turning to a final example it is instructive 
to consider the deep connection between ODE 
for Feynman diagrams, IBP identities and an- 
alytical properties of the diagrams. It can be 
seen as follows: for a given set of momenta we 
consider the transformation Pi = Tij(z)pj, sub- 
ject to 2 = E P = 0. Consider a general- 
ized, scalar, box (arbitrary powers in propaga- 
tors); we will also need the IBP equations for 
D (l, 1, 1, 1) and will define D (l) = D (l, 2, 1, 1) 
till D (4) = D (2, 1, 1, 1). Again, we can use IBP 
to get 

R4- tij D (j) = S i4 D (1,1,1,1) + A D (i), (13) 

where ADo(i) contains only 3 -point functions. 
Introduce the Caley determinant C4; it follows 
that P4 U = 2 M4, where {/ is unimodular (a sim- 
ilar relation holds for arbitrary N), i.e. 



dcti? 4 = 16C 4 , 



(14) 



so that the differential equation for the trans- 
formed box is 



d X. 
— D (1,1, 1,1) = — D (1,1, 1,1) 
dz C4 



Y. 



A straightforward calculation shows that 
1 d 



X 



2 dz 



C 4 



(15) 



for all values of {p}, {m} and for an arbitrary 
transformation T. Eq. (fT4]) holds for all TV, i.e. 
det R N = 2 N C N , where one should remember 
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that in four dimensions C N = for N > 6. The 
homogeneous term has the general stucture 

-i dC N 



f/2 



+ (AT - n) B^ 1 



dz 



where Y N _x is a combination oi N — 1 -point in- 
tegrals. 

5. An explicit example 

Our last case in point is given by the ODE for 
H — > g(pi)g(jP2) decay amplitude. Here there is 
one form factor F D that can be written, without 
reduction, as F D — J2i Fit 



{n-2)F 3 



An f Ml -2m* 
2 J 9 [0][1][2] 



= -2A„ / d™ 9 



— A r , 



q-pi 

'[0][1][2] 



16 



rf™<7 



q-piq-P2 



(6-n)q 2 



(16) 



Suppose that M 2 < 4 m 2 : define the transforma- 
tion Pi — Tij pj with 



T 



z 1 — z 
1 — z z 



Then S -» A/ 2 C/G with C = r 2 + /tif (1 + Ar) 
andG = -\M\ (l + 4r), being r = z(z-l) and 
u, 2 M 2 H = m 2 . The solution for AT is given by 



co < r AT < 




(17) 



The corresponding system of ODE will be written 
in terms of F\^ and F D giving 



dr 



F i =X ij F j +Y j , i,j = 1,2, D, 



(18) 



where X and Y are obtained by using IBP tech- 
niques, e.g. 

2 



X D \ — X DL 



1 + 4r 



X r 



0. 



(19) 



etc, with Y given by generalized two-point func- 
tions. Boundary conditions at AT are 

[M 2 

F D at (l + 6r AT ) 

- m 2 (l + 4r AT )] C$°*(z), 
F 1 Z T l -{M 2 H -2m 2 ) C«**(z) 

F 2 at M 2 z AT Ct S (z). (20) 

More general transformations, not singular for 
any z S R, exist but will not discussed here. The 
pure-scalar term becomes 



C Q (r)=Cr s (r)+cl e Hr) 



^ sing 



(r) In 



M 2 



Co S (r) 



1 + Ar 
with boundary conditions 



(21) 



<ru,-r, = gi/?(r AT ) /3 2 (r) = l-4^ 



H 



while the regular part is computed numerically 
(boundary condition for the regular part will not 
be reported here). 

The general strategy, e.g. for processes with 
A" = 4, is as follows: define 

D no ... nM A„ (d ^S V-, (g -T (22) 



[0] 



[3] 



which satisfy a recurrence relation (IBP) 

Bno.-nsW = ^ 1 rfn ...n 3 (i)+< ... Tl3 W, (23) 

then find the minimal set of linear combinations 
F = cD such that Amp — J2 F with {F} closed 
under d/dz. 

6. Extension to multi-loop 

Although we shall not discuss higher loops in 
details here, we present one simple example: the 
equal mass two- loop sunset S [ [T7] ; with scaled 
masses m = 1 and p 2 = x we perform the 
transformation i^zi 

d 2 

x z (xz + 1) (x z + 9) S(x, z) = 
dz 1 

P(x, z) — S(x, z) + Q(x, z) S(x, z) + R(x, z) 



dz 
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The AT solution is z AT = — x~ 1 (note that here 
AT = pseudo-threshold). For different masses we 
map 

Z % AT 1 ^ / \ 

nii -> Mi = — m l + — m, (24) 

1 Zat t ^ AT 

and use the previous calculation of AT. 

7. Conclusions 

In conclusions we have presented a proposal 
for solving the problem of reducing Feynman di- 
agrams which is based on a single variable defor- 
mation of the amplitude. 
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